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Abstract 



We present a new algorithm for the reduction of one-loop tensor Feynman integrals within the 
framework of the XLOOPS project, covering both mathematical and programming aspects. The 
new algorithm supplies a clean way to reduce the one- loop one-, two- and three-point Feynman 
integrals with arbitrary tensor rank and powers of the propagators to a basis of simple integrals. 
We also present a new method of coding XLOOPS in C++ using the GiNaC library. 
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1 Introduction 



The purpose of this paper is twofold. First, we introduce a new method to code XLOOPS in C++ 
with the help of the GiNaC library [1]. This also serves as an introduction to GiNaC for symbolic 
computations. 

Second, we introduce a new procedure for tensor reduction which allows us to manipulate one- 
loop one-, two- and three-point Feynman integrals with arbitrary tensor rank and arbitrary powers 
of propagators. Working in orthogonal and parallel space of momentum configuration [2], this proce- 
dure allows one to reduce the one-loop tensor integrals to a basis of simpler integrals without using 
the Passarino-Veltman procedure. This procedure has been implemented in the new version of the 
XLOOPS program, called XLOOPS-GiNaC. 

Unlike the original XLOOPS, XLOOPS-GiNaC currently only handles the one-loop case (and as 
such competes with existing programs for one-loop calculations such as those presented in [7]) but it 
has been designed with being a prerequisite for two-loop problems in mind, and the authors' goal is 
to make XLOOPS-GiNaC a powerful tool for one- and two-loop analytical Feynman integrations. 

In sections 2 and 3 we introduce briefly the motivation to rewrite XLOOPS in C++, and the 
GiNaC library. In section 4 the procedure used to reduce one-loop two- and three-point integrals is 
presented. In two appendices simple examples of C++ programs using the GiNaC library and the 
XLOOPS-GiNaC interface for one-loop integrals are presented. 

2 XLOOPS-GiNaC: The Motivation 

In the past, the XLOOPS package [6] has been developed in a heterogenous environment: The core 
routines for transforming Feynman graphs into the basic integrals and for analytic integration are 
implemented in the language of the Maple [8] computer algebra system, the graphical user interface 
is written in Tcl/Tk [9], and the numerical integration is done by C programs that are generated and 
compiled at run-time by XLOOPS. 

This way of implementation has a couple of drawbacks: 

• While Maple and other computer algebra systems provide sophisticated mathematical capabil- 
ities, it is not suited as an environment for developing large applications such as XLOOPS, as 
it was not primarily intended as a programming language and only offers limited support for 
modern software engineering. For example, the only structured data type in Maple is the list, 
the distinction between local and global variables is not consistently enforced by all language 
constructs, and tools such as debuggers are very rudimentary. 

• Different versions of Maple are in places incompatible with respect to the language. This leads 
to XLOOPS having to provide two versions of all program parts written in Maple one for Maple 
V Release 1 and one for Release 3. With Release 4 and later releases XLOOPS cannot currently 
be run. Not only does the parallel development of multiple program versions require higher 
maintenance effort, but it is also impractical for the user to require the installation of a specific 
version of Maple just to run XLOOPS. 
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• Maple is a commercial system, which prevents XLOOPS from getting a wide reach of distribution, 
especially in academic institutions. 

• The communication between different program parts is difficult because Maple has insufficient 
support for embedding programs written in other languages, it requires higher efforts for the 
conversion of data and, in some cases, for the multiple storing of redundant information (for 
example, the information about graph topologies has to be duplicated for the calculation and 
for the graphical interface). 

• In general, it is more difficult to maintain and develop a heterogenous program package. For 
the programmer it is necessary to become acquainted with three different environments (Maple, 
Tcl/Tk, and C). 

The listed drawbacks and the observation of bugs related to the internal structure of XLOOPS 
led to the conclusion that the program had reached a state in which further development was almost 
impossible (approx. 15000 lines of badly documented source code, in some places very obscure pro- 
gramming techniques like identifiers whose meaning depends on capitalization). Thus, the decision 
was made to rewrite XLOOPS from scratch, putting it on more solid grounds. 

In particular, it was decided that the new version of XLOOPS should be written in one uniform 
programming language. As the essential part of XLOOPS are the analytical calculations, traditional 
computer algebra systems (Maple, Mathematica [10], Reduce [11] and MuPAD [12]) were envisaged at 
first. But all these share more or less the same deficiencies as Maple, especially the low suitability of the 
built-in language for the development of large systems. It was therefore decided to use an established 
programming language as the foundation and extend it by the required algebraic capabilities. During 
the development of XLOOPS it was noted that only a small subset of the mathematical functions 
provided by Maple are actually needed. These are: 

• complex arithmetic with arbitrary precision; 

• simple manipulation of symbolic expressions, like expansion, collection, substitution of variables; 

• simplification of rational functions; 

• symbolic differentiation; 

• special functions like polylogarithms and the Gamma function; 

• Laurent series expansion of expressions containing these functions; 

• solving systems of linear equations for the analysis of the tensor structure of the integrals; 

• handling expressions containing elements of some particular non-commutative algebras such as 
Clifford and SU(3) Lie algebras; 

• numeric integration. 

Features that are particularly not required are: 
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• symbolic integration since only few master integrals need to be really integrated and they will 
be done by hand; 

• calculation of limits; 

• treatment of domains and assumptions about the range of values of variables. 
Due to the following reasons, C++ has been chosen as the programming language: 

• C++ is officially standardized [13], so fewer complications are expected from the future develop- 
ment of the language. 

• C++ allows to write down symbolic expressions in their natural mathematical notation by means 
of operator overloading (e.g. 4*a+b instead of something like add(mul(4,a) ,b)). 

• C++ is available for virtually all computer platforms. In particular, there are free compilers for 
the Unix systems predominant in the academic area. 

• There is a large assortment of development tools available like powerful source-level debuggers 
and systems for version control and documentation. 

• There is also a large number of existing libraries, especially for arbitrary precision arithmetics. 

• As a compiled language, C++ is also suitable for numeric integration. 

Based on C++ we developed a system called "GiNaC" 1 that is primarily aimed at the re-imple- 
mentation of XLOOPS, but is also suited for the development of other systems that integrate algebraic 
and numeric calculations with user interfaces for methods of data acquisition. 

3 GiNaC: A New Programming Environment For XLOOPS 

GiNaC [1] is a C++ library for handling symbolic mathematical expressions 2 . Some of the features of 
the library are 

• complex arithmetic with arbitrary precision, based on the CLN library by Bruno Haible [14]; 

• manipulation of symbolic expressions; 

• normalization of rational functions; 

• matrices and systems of linear equations; 

• numerous special functions (trigonometric and hyperbolic functions, exponential functions, log- 
arithms, Gamma und polygamma functions); 

• symbolic differentiation; 

1 an acronym for "GiNaC is not a Computer Algebra System" 
2 http : / /www . ginac . de 
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numeric 


real and complex numbers (integers like 42, exact frac- 




tions like | and floating point numbers like 7.319), 




based on CLN 


symbol 


algebraic symbols 


constant 


constants like tt which are treated similarly to symbols 




but also have a predefined numerical value 

n 


add 


sums of the form > axi, for a € C and arbitrary 




i=l 

non-numeric expressions Xi 

n 


mul 


product of the form X{ C% , for a € C and arbitrary 

: -1 




% — 1 

non-numeric expressions Xi 


power 


arbitrary expressions of the form x y 


pseries 


compact representation of Laurent and Taylor series 




(only contains the series coefficients, the expansion 




variable, and the expansion point) 


function 


mathematical functions like sin() und cosQ, where the 




individual functions are not implemented as subclasses 




of function but are distinguished by a function index 


1st 


lists of expressions 


matrix 


matrices 


relational 


equations and unequations 



Table 1: The most important GiNaC classes 

• series expansion of functions (Taylor and Laurent series); 

• Clifford and SU(3) color algebras (this is a recently added feature in GiNaC, making it especially 
suited for applications in particle physics); 

• it is Open Source, licensed under the GNU General Public License. This means that it is not 
only freely available but also free in the sense that users have access to the sources and are 
allowed to modify or extend the library and to redistribute the modified version. This is in 
sharp contrast to most other CAS that place heavy restrictions on their legal use. 

GiNaC is designed in an object-oriented fashion. The central class of GiNaC is the class ex that 
stores a symbolic expression. Strictly speaking, ex only represents a "smart" pointer to the real 
expression which is stored as a tree whose nodes are objects subclassed from the abstract base class 
basic. The operators +, -, * und / are overloaded to simplify the creation of expressions in the 
program code. 

Table 1 gives an overview of the most important subclasses of basic. The classes can be categorized 
into the "atomic" classes numeric, symbol and constant which are at the leaves of an expression tree, 
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and the container classes (all others) which themselves contain expressions. The representation of 
sums and products with numeric coefficients and powers was chosen for reasons of efficiency (see [15] 
and [16]). 

A detailed description of the internal functionality of the basic classes and methods of GiNaC is 
given in [1] and [16]. 



4 One-Loop One-, Two-, And Three-Point Integrals 

In this section we present our algorithm for tensor reduction, which was introduced by Collins [4] and 
further developed by Kreimer [2, 3]. A similar approach was used in a subsequent paper by Ghinculov 
and Yao [5] . The advantages of this approach compared to other procedures like the Passarino-Veltman 
procedure, especially with regard to future two-loop applications, are outlined in [16]. 

To regularize UV-divergences, we are working in D = 4 — 2e dimensional space-time. The main 
idea is to split the space of integration in a parallel and an orthogonal space. We define the parallel 
space to be the linear span of the n-external momenta qi (i = 1, n) involved in the integrand. This 
parallel space has a finite dimension J < D. The remaining D — J dimensions span an orthogonal 
complement of the parallel space called the orthogonal space [2]. 

Once an explicit configuration of external momenta is chosen, the dimension of the parallel space 
J is known, and the scalar products are written explicitly in terms of the components of the external 
momenta 

;2 _ ;2 ;2 ;2 j2 

I — t Q — l 1 — ■ ■ ■ — — l ± , 

I ■ Qi = loQi - hQh — ■ ■ ■ — Ij-iqij^. 
A general one-loop integral can be written as 

£} — J oo oo oo 

J d D lF(l 2 ,{l-q t }) = f ^-^j^ j dl --- j dlj-! Jdl ± l I l- J - 1 F(l ± ,l ,...,lj- 1 ). (2) 

— oo — oo 

In the usual Passarino-Veltman approach, to keep an explicitly covariant form, one expands tensor 
integrals of the type 



in a basis of Lorentz tensors constructed from the metric tensor g^ v and the external momenta qi^. 

(n) 

Since Tf ll ',^ N is symmetric, one can re-group indices belonging to the parallel space and to the or- 
thogonal space together 

ji(n) _ yr(re) 

fll—HN flj...fl D 

PO Pl PJ-1 
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where po,Pi, ■ ■ • denote the numbers of indices for the 0-, 1-, ... components. The indices fij, . . . /id 
in the orthogonal space have to result in a symmetric combination of metric tensors in the orthogonal 
space g^ v 



D-J 



T (n) = Y 1 ) 2 ( ± a ± \ T {n) (5) 

■ i /*l.../ij V J£ \ y ^.TfJ.J+l ■ ■ ' y ^D-lHD J symm (P0-PJ-1PX)" V ' 

with N = po + ... +p±. The normalization can be derived by looking at contractions with g^ v 

and observing that 

ai^y^ 1 = d - j. (6) 

One obtains 



k= n ip-j+2%). (7) 

i=J 

The coefficients needed for the calculation of a specific component of a general one-loop tensor integral 
therefore have the form 

/jPO ]PJ-iiP± 
( f> l L ■■■ L J-1 l ± /gx 

rr™ /Y i \2 2 , ■ Y 
11=1 + H) -mf + ip) 

(n) 

In the next sections, we present an algorithm to calculate TY p ^ i n for the one- loop two- and 

(n) 

three-point tensor integrals. In the rest of this paper, we call p ^ tensor integrals. 
4.1 One-Loop Two-Point Tensor Integrals 

In this section, the algorithm for an automatic calculation of one-loop two-point tensor integrals is 
presented. We first consider the case where q 2 is timelike, q 2 > 0. Then one can choose a reference 




l 

Figure 1: Notation for two-point functions 

frame where q^ = (qo, 0, 0, 0). The general integral for a tensor Feynman diagram shown in Fig.l has 
the form 

^(^ = ^1^ (9) 
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with 



Pi = (k + qo) 2 ~ l 2 i_ ~ m\ + ip, 
P2 = ll-l\-ml+ ip, 



(10) 



where In and l± span the parallel and orthogonal subspaces respectively and the integral vanishes 
unless j is even. In the spacelike and lightlike cases where q 2 < or q 2 = one can choose a reference 
frame where = (0, qoi,0, 0) or q^ = (501,901)0,0) respectively and the integral space can be split 
into a two-dimensional parallel and aD-2 dimensional orthogonal subspace. We will consider these 
cases later in section 4.2.3. A genuine one-loop integral has t\ = t 2 = 1, but the more general case is 
needed in the case of the reduction of integrals with more than one loop. 

The strategy now is to express I^ 2 J^ (q 2 ) as a polynomial of simpler integrals. It turns out that 
the usual scalar one- and two-point integrals An and E>n (in Passarino-Veltman notation [18,19]) are 
sufficient. This expansion is always possible except for some special cases that we will consider later. 

Firstly consider the general case where go 7^ 0- We express the numerator of the integral in Eq.(9) 
as a function of Pi and P 2 . From Eq.(10) we get: 



In -> l (P 1 ,P 2 ,q ,m 1 ,m 2 ) 

l\ -> l 2 L (P 1 ,P 2 ,qn,m 1 ,m 2 ) 
Ci 

C 2 = -ml + ip. 



2<?o 



(Pi 



Ci), 



In — P2 + c 2 , 



(11) 



2 2,2 
Qo ~ m l + m 2> 



Inserting Eq.(ll) into Eq.(9) and expanding the numerator of the integrand, one obtains 



At\<i 2 ) = Y:c nm jd D ipr H p, 

n,m 



m—t2 
2 



(12) 



with C nm being simple functions of qn, mi, m 2 , and 

r 



jd D ip[ 



n—t\ T}in—t2 



pt2-m 
r 2 



pn—ti 

d D r i 



pt 2 -m 
r 2 



/T>m—t2 
pti-n 



if n — t\ > and m — t 2 > 0, 



if n — t\ < and m — t 2 < 0, 



if n — t\ > and m — t 2 < 0, 



if n — t\ < and m — t 2 > 0. 



(13) 



We see that the second case actually corresponds to a scalar two-point function. For the last two 
cases, from Eq.(10) one can insert 



Pi = Pi(P2, In, qo) = 2l q + P 2 + Ci, 



(14) 



8 



or 



P 2 = Pi{Pi, h, 9o) =Pi~ 2l q ~ Ci (15) 

and expand the numerator of the integrands in Eq.(13). Then Eq.(lO) can be completely reduced to 
the one-point functions 



[ d Dj 

(P — m 2 + ipY 



4.1.1 The case q = 

If qo = 0, Eq.(13) must be rewritten using 



that is actually the one-loop one-point function. 

If mi 7^ one performs partial fraction decomposition and finds 

4tV)= fdPi lhf± 



tlt2 w 1 J P 1 (m 2 1 ) t i p 2 (m 2 2 y* 



1 ( 1 r- p n 



(ti - 1)! d(mf)*!- 1 Umf-mD^y Pi(mf 



(16) 



(17) 



-Pi = l — l± — m 1 + ip, 
P2 = Iq - l± - m 2 , + ip. 

If mi = m2 then Pi = P2 and 1^,^ {q 2 ) has the simple form 

e ? V)=/^j^ (is) 



(19) 



(t 2 - 1)! dCml)**" 1 ^(ml-mf)*! 7 P 2 (m|) J ' 
which is a combination of one-loop one-point functions. 

4.2 One-Loop Three-Point Tensor Functions 

The notation we use for one-loop three-point functions is shown in Fig. 2. We are working in the frame 
of reference where the external momentum configuration is q\ = {qw, 0,0,0), qi^ = (q20, Q21, 0, 0). 
The parallel space is now two-dimensional and the general form of the one-loop three-point tensor 
function is 

/■(3)«fc_ />; l l l l ± (on) 
1 h t 2 t 3 - J a L p ti p t 2 p t 3 V ZU J 
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Figure 2: Notation for three-point functions 



with 



-Pi = l l ~ l i ~ l ± + 2 ^0 9io + 9io -m\+ ip, 

,2 ,2 



-F*2 = ^0 ~ ~ l ± + 2 ^0 920 - 2 921 + 920 ~ 921 ~ m 2 + ¥, ( 21 ) 

^3 = ll ~ l\ ~ l\ ~ m 2 3 + ip 
where {Iq, h} and l± span the parallel and orthogonal subspaces, respectively. 

4.2.1 The general case 

Firstly we consider the case where q±o ^ and q2i ^ 0. Then we are always able to express Iq, Zi, l± 
in terms of P±, P2, P3: 

k = 7^-(Pl-P3-C 00 ), 

2 9io 

h = 2qw q 21 [ Q2 ° ( Pl ~ P 3 ~ («10 - m l + m D) + <?10 ( P 3 ~ P2) + 910 {m\ -m 2 2 + ql - q^)] 
= 7^— [ciolo + (P3 ~ P2) + cu] , 

^921 

l\ = Iq - l\ — P3 - C20 
with 

coo = 9io ~m\ + ml, 

cio = 2 920, 

2 2,2 2 

en = m 3 -m 2 + q 20 - q 21 , 
c 2 o = ml - ip. 

Again, as in the case of two-points functions, one substitutes Eq.(22) into Eq.(20) and obtains a 
combination of scalar three-point functions 

4tf = E c ™ d D ip[ m - H) pt t2) pt tz) m 

m,n,r 
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with C mnr being simple functions of masses and components of external momenta. More explicitly, 
the integrand on the right hand side of Eq.(24) contains terms like 

' and WT^T (25) 



pni r>W2 p«3 ' r> n k ' r> n j r> n k 

^1 r 2 ^3 

with i ^ j ^ k and positive rn. Other possible combinations lead to a vanishing integral. 

The first group of terms can be obtained from derivatives of scalar three-point functions. For 
the last two cases, using Eq.(21) one can expand the numerator in terms of propagators P n . in the 
denominator. This expansion is always possible and reduces the second group of terms to one-point 
functions. Similarly the third group of terms can be reduced to one-loop two-point functions. Note 
that in this step we meet two kinds of one-loop two-point functions. The first kind is the one-loop 
two-point function with one parallel dimension that was already found in the previous section. The 
second kind is the one-loop two-point function with two parallel dimensions of the internal momentum 
I that is not trivial and will be given in the section 4.2.3 as a separate case. 

4.2.2 The case q w q 2 \ = 

In this case the expansions in Eq.(22) cannot be used. However, Eq.(20) can be reduced by partial 
fraction decomposition. First consider the case qio = with arbitrary values of q 2 i- 

4.2.2.1 The case q w = 

Eq.(21) simplifies to 



(26) 



If mi ^ m3, partial fraction decomposition leads to a separation into terms which contain only 
two propagators 



Pi 


- I 2 


-If 


-ll 


-m 1 + ip, 


P2 


— '0 


-If 


-ll 


+ 2Z <?20 - 2hq 2 i + g| - q\ x -m 2 2 + ip. 


P3 


- 1 2 


-If 


-ll 


— m 2 + ip. 



T (3)ijk TT ( l_ d'l 1 

\}\{t f -l)\d{m 2 ff-' 




If m\ = 7713, two propagators are equal and it is sufficient to calculate 

T (3)ijk _ f,D, ^0 A l ± / 9 o\ 

tlt2 ' 3_ r Pl^{m 2 )P?{m 2 ) [ ] 

The three-point integrals are then reduced to combinations of two-point integrals with two parallel 
dimensions. We will treat this class of two-point functions in section 4.2.3. 
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4.2.2.2 The case q w ^ and g 2 i = 

In this special case, the integrals collapse from two parallel dimensions to one and Eq.(21) reads 

Pi = (k + qiof -l'l-m\+ ip, 

P2 = (Zo + <?2o) 2 - l'l ~ m\ + ip, (29) 
P 3 = lo-l'l-m 2 3 + ip 

where the components l\ and l± can be combined into l'j_ to form a new D — 1 dimensional orthogonal 
subspace with 

If = l\ + l\. (30) 



l£li+ k = d D l ° l j f . (31) 

f l f 2 *3 / p*l p*2 p*3 



Then the integral can be reduced to a simpler integral with only one parallel dimension 

M 

pi'p^Ps 

Again, one can use the same procedure as in the general case 

l -> lo(Pi,P3,m,qio) = ~ — (-Pi - P3 + m\ - ml - qf ), 

•%o (32) 

l± l±(Pi,P3,mi,qio) = Iq - P3 - ml + ip 

that reduces Eq.(31) to a form similar to Eq.(24). 

4.2.3 The two-point integral with two parallel dimensions J%* k 

In the preceding section we have shown that general one-loop three-point tensor functions can be 
reduced to the usual scalar integrals and one new two-point function corresponding to a tensor com- 
ponent in a two-dimensional parallel space. Explicitly, this integral reads 

jijk _ I ' jD, ^0 K f± 

1 [{lo + qio) 2 -ll-ll-m\ + ip]*i [(/ + qio? ~ (h + ?2i) 2 -l 2 ± -m* + ip}** 

•D, (k-qioY Ijlj 

[l 2 -l 2 -ll-m\ + ip]ti [(l + Qo) 2 - (h + Qi?-ll-ml + ip}** 

with Qo = q2o — qio, Qi = 921; = <Z2 M — If Qi = 0> this integral reduces to the one-loop 
two-point function in one-dimensional parallel space as found in the previous section. If, on the other 
hand, Q\ 7^ 0, one can always find a Lorentz boost which transforms into a reference frame where 
the transformed 4-momentum Q'^ has either only one non-zero component (Q' or Q\ if Q is timelike 
or spacelike) or where Q' = Q[ if Q is lightlike. The loop momentum has to be boosted accordingly 
which, however, modifies only the numerator of the integrand in Eq.(33). Explicitly, consider the 
boost 

l o\ _ ( 7 10\ fl'o 



t 1 , , 1 1 ,/ 1 (34) 



Then the three sub-cases are treated as follows. 
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4.2.3.1 The timelike case Ql - Q\ > 

In this case, under the transformation in Eq.(34) the integral will be reduced to one- loop two-point 
functions with a one-dimensional parallel space as found in the previous section: 



juk = i d D lt b(i'o + Pii)-qi] 1 hWi'o + i'iW i± 

' [Itf-lf-l'jj-ml + ip]* 1 [(l> + Py-l[ 2 -l>f-mj + ip] t2 



with P = y/Q 2 - Ql 7 = Q /P and (3 = Qi/Q . 

4.2.3.2 The spacelike case, Ql - Q\ < 

In this case the boost with P = \J Q\ — Q 2 ,, 7 = Qi/P and (3 = Q0/Q1 transforms to a reference frame 
in which the integral reads 

jijk = f d D lf [7^0 + W-gir h(Pi' + i[W i't f36) 

2 J [Itf-lf-l'l-ml + ip]* 1 [l^-{l> l + PY-l'l-ml + i P ] t2 

The components 1' and l' ± can be combined to form a new D — 1 dimensional orthogonal subspace 
while l\ spans the parallel subspace. Using the same procedure as in the previous sections one can 
reduce the integral completely to the scalar one- and two-point functions. 

4.2.3.3 The lightlike case, Q 2 - Qj = 

In this case the transformation in Eq.(34) is singular and the integral J^ k becomes 

jijk _ f, D] (*0 ~ gl)' l l l l /o 7 x 

J 2 — u 1 u2 72 r2 __2 , - if, r/7 , /-> \9. n , ^> \2 i2 __2 , -Uo' V° ' ) 



(38) 



[ll -il-ii-m\ + ip]*i [(/o + Qq? - (h + Qo) 2 - ll - ml + ip\* 
By solving the system of equations 

p i = l l ~ % ~ l± ~ ml + ip, 
P2 = (lo + Qo) 2 ~ (h + Qo) 2 -l\-ml + ip, 
one obtains 

l\ = ll~l\- P\-m\ + ip, 

, P l -P 2 + ml-m{ , (39) 

h ~ Wo + 

Inserting Eq.(39) into Eq.(37) and expanding the numerator of the integrand, the integral will be 
reduced to scalar one-point functions and a simpler tensor integral 



Ji = ld D l ^ i If;. 

1 ^ [ll -ll-ll-m\ + ip]«i [(l + Qo) 2 - (h + Qo) 2 - l 2 ± - ml + ipf* ' 

The explicit calculation of this integral is given in [17]. 

This completes the description of our algorithm for tensor reduction. We did not reproduce explicit 
expressions for the basic scalar integrals in this paper since these can be found in the literature [2, 19]. 
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5 Conclusion 



Due to the limitations of Maple and the internal structure of XLOOPS we decided to rewrite XLOOPS 
from scratch, based on GiNaC, a C++ library to replace Maple as an algebraic programming envi- 
ronment [1]. An efficient algorithm for one-loop one-, two- and three-point tensor reduction was also 
successfully implemented. At this stage of the project, a package for calculating one-loop one-, two- 
and three-point tensor integrals is available and can be downloaded from http://wwwthep.physik. 
uni-mainz.de/~xloops. Like GiNaC, it is distributed under the terms of the GNU General Public 
License. 

As the next step, we plan to rewrite the module for two-loop one-, two- and three-point integrals 
and to completely recode XLOOPS in C++ using the GiNaC library, providing a package for doing 
particle physics in a homogenous C++ environment. 
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Appendix A 

In this section we introduce the definitions for the one-loop integral functions in XLOOPS-GiNaC. 
The XLOOPS-GiNaC package provides the three GiNaC functions OneLooplPt () , 0neLoop2Pt () , and 
0neLoop3Pt () , which can be used in algebraic expressions like any other predefined GiNaC function: 

• The one-point function 




• The two-point function 





[(lo + q) 2 -ll-m\ + ip]*i [ll -l\-ml + ip]* " 



• The three-point function 



DneLoop3Pt(i, j , k, q 10 , q 2 o, q2i, mi, m 2 , m 3 , t l5 1 2 , t 3 , 



r(3) ijk 
<1 <2 *3 
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with 



Pi 


- J 2 


-q-il 


+ 2/ qio + qf -m\ + ip, 


P2 


- I 2 


-H-H 


+ 2l q 2 o-2 h q 2 i + q% - q\\ - ™| + ip. 


P3 


- 1 2 




— m 2 + ip. 



As with any GiNaC function, the arguments and return values of the above functions are objects 
of type ex so the return values as well as input parameters can be any symbolic or numeric expression. 

In order to illustrate the output of XLOOPS-GiNaC, we give one example program that calculates 
and prints out both analytical and numerical results of the UV-divergent and the finite terms of the 
integral 0neLoop2Pt(l, 0, q, mi, m 2 , 1, 1, p) 

1 #include <iostream> 

2 #include <ginac/ginac .h> 

3 #include <xloops/xloops .h> 

4 using namespace GiNaC; 

5 using namespace xloops; 

6 

7 int mainO 

8 { 

9 symbol q("q"), ml ("ml"), m2("m2"), epsC'eps"), rhoO'rho"); 

10 ex a = DneLoop2Pt(l, 0, q, ml, m2, 1, 1, rho) ; 

11 a = a.series(eps == 0, 4); 

12 ex al = a.coeff(eps, -l).subs(rho == 0); 

13 ex a2 = a.coeff(eps, 0).subs(rho == 0); 

14 cout << "Order eps~-l is " << endl << al.normalO << endl; 

15 cout << "Order eps~0 is " << endl << a2.normal() << endl; 

16 ex bl = a. subs (rho==0) .subs (ml==80) .subs (m2==80) .subs (q==100) .evalfO ; 

17 cout << "Numerical value up to order eps~2 is " << endl << bl << endl; 
is return 0; 

19 } 

The output of this program reads 

1 Order eps~-l is 

2 -l/2*I*Pi~2*q 

3 

4 Order eps~0 is 

5 (-l/2*q"(-l)*m2"2-l/2*q+l/2*q"(-l)*ml~2)*(I*Pi~2*q~(-l)*((l/2*q-l/2*q~(-l)*(ml 

6 -2-m2-2))*R2exl(-m2-2,-(l/2*q-l/2*q-(-l)*(ml-2-m2-2))-2)+R2exl(-ml-2,-(l/2*q+l 

7 /2*q"(-l)*(ml"2-m2"2))~2)*(l/2*q+l/2*q"(-l)*(ml"2-m2"2)))+q*(-I*Pi~2*log(Pi)*q 
s "(-l)+Pi"(3/2)*(I*sqrt(Pi)*q"(-l)*(-2*log(2)-Euler)+2*I*log(2)*sqrt(Pi)*q"(-l) 
9 +2*I*sqrt (Pi) *q~ (-1) ) ) -Pi~3)+l/2*q~ (-1) *(-I*Pi~2*m2~2*log(m2~2) -I*Pi~2*log(Pi) 

10 *m2"2+Pi"2*(I-I*Euler)*m2~2)-l/2*q" (-l)*(Pi"2*(-I*ml~2*log(ml~2)-I*Euler*ml~2+ 

11 I*ml~2)-I*Pi~2*log(Pi)*ml~2) 
12 

13 Numerical value up to order eps~2 is 

14 (-493 . 48022005446793098*1) *eps~ (-1) + (5019 . 9161138633880865*1) + (-4 . 6074255521 

15 943996428E-15-25944. 085010687200793*1) *eps+0rder (eps~2) 
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Appendix B 



In order to illustrate the use of XLOOPS-GiNaC for calculating one-loop Feynman diagrams, we give 
one example program that is actually part of the automated regression tests of the package. It checks 
that the longitudinal part of the vacuum polarization in QED vanishes on the one-loop level as required 
by gauge invar iance, i.e. 

<Z M IPV) = 0, (41) 

where 

The general tensor structure of IP" is 

q 2 

with functions A and B that, because of (41), satisfy 

A + B = 0. 

This expression is obtained by contracting IP" with < ^ L : 



wSe*L = A + B 

r 

o2 



= - e l (d D l Tr 4d + 4 + ™M + m) 
q 2 J {{l + q? -m 2 )(/ 2 -m 2 )' 



(42) 



The following program verifies that this expression vanishes for the first three orders of the regular- 
ization parameter: 

1 #include <ginac/ginac .h> 

2 #include <xloops/xloops .h> 

3 using namespace GiNaC; 

4 using namespace xloops; 

5 

6 int main (void) 

r { 

s symbol D("D") ; 

9 symbol 1("1"), q("q"), m("m"), e("e"); 

in symbol 10("10"), lorthC'lorth") , epsC'eps"), rhoC'rho"); 

n 

12 scalar_products sp; 

13 sp.add(l, 1, pow(10, 2) - pow(lorth, 2)); 

14 sp.add(q, q, pow(q, 2)); 

15 sp.add(l, q, 10*q) ; 

16 

17 ex I = -pow(e, 2) / pow(q, 2) 
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* dirac_slash(q, D) 

* (dirac_slash(l , D) + dirac_slash(q, D) + m * dirac_DNE 

* dirac_slash(q, D) 

* (dirac_slash(l, D) + m * dirac_DNE() ) ; 
I = dirac_trace (I) . simplif y_indexed(sp) ; 

ex a; 

for (int i=0; i<3; i++) 

for (int j=0; j<3; j++) { 

ex c = I.coeff(10, i) . coef f (lorth, j); 

a += c * DneLoop2Pt(i, j, q, m, m, 1, 1, rho) ; 

} 

a = a.series(eps == 0, 4); 

ex al = a.coeff(eps, -l).subs(rho == 0); 
ex a2 = a.coeff(eps, 0).subs(rho == 0); 
ex a3 = a.coeff(eps, l).subs(rho == 0); 

cout « "Order eps~-l is " « al.expandO « endl; 
cout « "Order eps~0 is " « a2.expand() « endl; 
cout « "Order eps~l is " « a3.expand() « endl; 

return 0; 
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Line Explanation 



1-4 Include the header files for the GiNaC and XLOOPS libraries. 

8-10 Declaration of all appearing symbols. These are the spacetime dimension D, the loop and 
external momenta I and q, the mass m, the parallel and orthogonal space loop components 
Iq and l±, the dimensional regularization parameter e, and the infinitesimal imaginary part 
of the propagator p. With GiNaC it is necessary to specify the name used for printing 
expressions because C++ does not provide the names of variables at run-time. 

12-15 The possible scalar products of the momenta are expressed in terms of q, Iq and l± and 
are registered in a scalar_products object which is later passed to s impl if y_ indexed () : 
l-l = Iq— q-q = q 2 ,l-q = loq. pow() is used for exponentiation because the C++ operator 
'has the wrong precedence in relation to the operator *. 

17-21 The numerator of the integrand of Eq.(42) is constructed using Clifford algebra objects in 
a nearly 1-to-l translation of the right-hand side of that equation. Note that GiNaC does 
not require the use of a special operator for non-commutative products here. 

22 The trace is taken and the resulting expression which contains metric tensors is sim- 
plified, inserting the scalar products defined above (the result for I is — 4e 2 (7g + l\ + 
qlo + m 2 )). Expressions are usually manipulated in the C++ oriented notation expres- 
sion, function (parameter 's) , but the functional notation function (expression , parameters) 
is also available, as shown. 

24-29 The integral is now expressed in terms of the basic 0neLoop2Pt() integral functions by 
assembling the coefficients of all powers of Iq and l±. 

30 To get the UV divergence as well as the finite part, one needs to take the series expansion 
of the integral at the pole e = 0. 

32-34 The coefficients of the series for the orders e^ 1 , e° and e 1 are extracted and the limit 
p — > is taken by calling the function subs(). 

36-38 The three coefficients are simplified with expand () and printed to the standard output 
stream. When run, this program will output all three coefficients as 0. 
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